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ABSTRACT 

It has been observed that Cepheids in the Magellanic C louds have lower masse s for the same 
luminosity than those in the Milky Way. The model, from [Neilson fc Lester! (|2008), of pulsation- 
driven mass loss for Cepheids is applied to theoretical models of Cepheids with metallicity consistent 
with the Milky Way and Large and Small Magellanic Clouds. The mass-loss model is analyzed using 
the metallicity correction of the Period-Luminosity relation to compare the ratio of mass loss of 
Cepheids with lower metallicity to that of Cepheids with solar metallicity. It is determined that mass 
loss may be larger for the lower metallicity Cepheids, counterintuitive to radiative driving estimates. 
Also the mass-loss rates of theoretical Cepheid models are found to be up to 5 x 1CP 9 for Galactic 
Cepheids, 5 x 10~ 8 for Large Magellanic Cloud Cepheids, and 2 x 10^ 7 M & /yr for Small Magellanic 
Cloud Cepheids. It is argued that mass loss increases as metallicity decreases for Cepheids with 
periods less than 20 days and that mass loss decreases for longer periods. Assuming dust forms in 
the wind of a Cepheid at some distance, the infrared excess of the models is computed, finding the 
infrared brightness is approximately a magnitude larger due to mass loss. The infrared magnitudes 
are compared to recently published Period-Luminosity relations as a test of our predictions. 
Subject headings: Cepheids — circumstellar matter — Magellanic Clouds — stars: mass loss 



1. INTRODUCTION 

£f~) ' The study of classical Cepheid variables has been 

0^ . a cornerstone of stellar and Galactic astrophysics 
' since the discovery of th e Period-Luminosity relation 
(|Leavitt fc Pickeri"nglll912D . Observations of light curves 
, and period provide information for testing stellar struc- 
■ ture and evolution while the Cepheid Period-Luminosity 

OO ' relation makes them ideal standard candles. It has been 
by understanding the physical structure and mechanism 
of pulsation that Cepheids have become powerful tools, 

. £^ \ but there are still phenomena being discovered. 

One phenomenon was the recent discovery of infrared 
excesses from circumstellar material surrounding Galac- 
C3 ; tic C epheids (JKervella et al.l I2006t IMerand et all I2006L 
l2007f) . It has been suggested that the cause of the in- 
frared excess is dust f orming in a wind from the Cepheids 
(jKervella et al.l 120061 ) . There is additional evidence for 
mass loss, albeit circumst antial and wide ranging. Analy- 
sis of IRAS observations (|McAlarv fc Welchlll98alDeasvl 
1988) infers mass-loss rates up to 10 6 M©/yr for some 
Cepheids, although radio observ ations have placed up - 
per l imits of lQ- 9 -lQ- 7 M^ Ivr (I Welch fc Duricl fl98l . 
Also iRodrigues fe Bohm-Vitensel (fl992) determined an 
upper limit of M < 7 x 10~ 10 M Q /yr for SU Mus and 
iBohm-Vitense fc Lovd (|1994l ) found M«5x lO" 5 M /yr 
for I Car. Furthermore the discovery of circumstel- 
lar material surrounding I Car, Polaris, 5 Cep , Y Oph 
(jKervella et alJ 120061 : IMerand et aTl I2006L 120071 ) as well 
as RS Pup and SU Cas being assoc iated with nebulae 
(|Havlenl Il972l: iTurner fc Evans! Il98l ) provide more evi- 
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dence for mass loss. 

Theoretical evidence fo r mass loss based on h ydro- 
dynamics is presen t ed by IWillson fc Bowenl (11986) and 
iBrunish fc Willsonl (fl989l ). In iNeilson fc Lester! (|2008l . 
hereafter Paper I) we developed an analytic method to 
calculate mass-loss rates for Cepheids in which the mass 
loss is driven by the combination of pulsation and ra- 
diation. The method was applied to observed Galac- 
tic Cepheids where the mass-loss rates are predicted to 
range from 10 -10 -10 -7 M G /yr. This implies that over 
the lifetime of Cepheid evolution, typically 10 to 30Myr, 
a significant amount of mass can be lost. 

Mass loss is interesting in Cepheids not only because 
it produces infrared excess and may affect the period 
of pulsation, but it has also been proposed as a so- 
lution to the Cepheid mass discrepancy, which is the 
difference between mass estimates of Cepheids by stel- 
lar evolution models and by pulsation models. Histor- 
ically, Hie mass determinations differed by as much as 
40% (ICoxl 1 19801) . The intro duction of the OPAL opaci- 
ties ()Rogers fc Iglcsias 1992) brought a partial resolution 
to the m ass discrepancy problem by reducin g it to or- 
der 10% (jMoskalik. Buchler, fc Maroml 11992ft. however 
the discrepancy is still an issue. iSebo fc Wool (EUl) 
found a mass discrepancy of about 20% for Cepheids 
with periods ranging from 2.7 to 30 days in the Large 
Magallanic Cloud (LMC) cl u ster N GC 1850; likewise 
iBeaulieu. Buchler. fc Kollathl |2001) found similar re- 
sults for Cepheids from the OGLE database for the LMC 
as well as a larger mass discrep ancy for Smal l Mag- 
allanic Cloud (SMC) Cepheids. iBrocato et all (|2004j ) 
found masses of Cepheids in the LMC cluster NGC 1866 
to be smaller than predicted by evolution by almost 30% 
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for short period Cepheids. The only dynamical mass es- 
timates for Cepheids are from observations of Galactic 
Cepheids in binary systems; these are all found to have 
masse s less than that predicted by evolutionary ca lcula- 
tions (jBohm-Vitense et alJll99l lEvans et al.ll2006[ ). 

The current status of the m ass discrepancy for G alac- 
tic Cepheids was presented bv lCaputo et a l. (2005), who 
give the discrepancy as a function of period. For short pe- 
riods, the mass discrepancy is approximate ly 10% - 20% 
and this decreases with increasing period. iKellerl (|2008f ) 
argues against this conclusion based on an analysis of 
the behavior of the mass-luminosity relation of Cepheids 
at large masses, and he proposes the mass discrepancy 
is constant with respect to mass, about 20%. The re- 
sults differ because of the stellar evolution models used; 
iCaputo et all (|200l use stellar evolution models that do 
not include a prescription for mass loss on the main se - 
quence or on the red giant branch while the Keller (2008) 
results use stellar evolution models that do include mass 
loss. Mass loss affects the evolution of the massive stars 
M > 10A/q and causes a non linear mass - luminosity 
relation for Cepheids, the vari ous mass-lu minosity rela- 
ti ons are shown in Figur e 2 of IKellerl (|2008l ). 

iKeller fc Woodl (12002 ■ explore the mass discrepancy 
of LMC Bump Cepheids using non-linear pulsation 
models and find a mass discrepancy of about 20%. 
IKeller fc Woodl (|2006f) determined that the mass discrep- 
ancy for SMC Cepheids is about 20% for Cepheids with 
mass 5-7M©, while the discrepancy for LMC Cepheids 
is about 17%, determined by matching MACHO and 
OGLE observations with non-linear pulsation models 
with varying mass and me tallicity. This obs ervation has 
been verified for the LMC (|Testa et al.ll2007l), and it also 
agrees qualitatively with the analysis of iCordier et al.l 

pool . 

While mass loss is one possible cause of the mass 
discrepancy, a second suggested source is convective 
core overs hooting in the main seq uence progenitors of 
Cepheids (jHuang fc Weigertl 1 19831 ). There is no direct 
evidence of this in Cepheids other than the mass dis- 
crepancy, coupled with the argum ent that mass loss is 
negligible over the lifetime of stars (|de Jager et al.|[i~988l ; 
iGirardi et alJ 12000). However there is evidence for con- 
vective core overshooting from stellar isochrones fit to 
stella r cluster observation (for e x ample iMucciarelli et al.| 
' 20071; IChiosi fc Vallenaril 120071 iRosvick fc Vandenbergl 
19981 ). ISchroder. Pols, fc Eggletonl (119971) used observa- 



tions o f C Aurigae systems and lRibas. Jordi. fc Gimenei 
(2000) used observations of detached binaries to con- 
strain convective core overshooting. Representing con- 
vective core overshooting as Acco — oiHp, where Hp is 
the pressure scale height, the free parameter a produces 
a mat ch to the observations for a ss 0.1-0.3. iRibas et al.l 
(2000) find larger convective core overshooting for the bi- 
naries V380 Cyg, and HV 2274. This evidence is circum- 
stantial because the parameter describing convective core 
overshoot is obtained by matching models that might ig- 
nore other physical processes. 

The purpose of this article is to investigate the effect 
of metallicity on mass loss in Cepheids using the ana- 
lytic model from Paper I. In particular, we wish to ex- 
plore the potential for Cepheid mass loss in the Large and 
Small Magellanic Clouds. Understanding Cepheid mass 
loss in the Clouds will provide insight into the source of 



the mass discrepancy. The goals of this work are to dif- 
ferentiate how mass loss and convective core overshoot 
would contribute to the mass discrepancy, as well as to 
predict observational consequences of mass loss, such as 
infrared excess. We start this exploration by deriving 
the metallicity dependence of the analytical mass-loss 
formulation from Paper I in Section 2. In Section 3, the 
analytical model is applied to Cepheids in the LMC and 
SMC to predict the potential behavior at lower metal- 
licity. Mass-loss rates are calculated in Section 4 for 
theoretical models of Galactic, LMC and SMC Cepheids 
and the resulting infrared luminosities are predicted in 
Section 5. The dependence of the mass discrepancy on 
pulsation-driven mass loss is explored in Section 6. 

2. THE METALLICITY DEPENDENCE OF THE 
ANALYTICAL MASS LOSS RELATION 

In Paper I, an analytic solution for the mass-loss rate 
as a function of pulsation pha se was derived based on 
an ext ension of the method of ICastor. Abbott, fc Kleir] 
(1975). This method solves the momentum equation 
which balances the force of gravity with the forces due to 
continuum and line radiative driving and contributions 
due to pulsation and shocks generated by the pulsation. 
An approximate solution to the momentum equation is 
derived in Paper I for the mass-loss rate at some phase 
of pulsation, 
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The mass-loss rate is an implicit function of the pulsa- 
tion phase and the integral of Equation [T] with respect 
to phase over the period is the phase-averaged mass-loss 
rate. The variables k and a are force multiplier param- 
eters that depend on the gravity and effective temper- 
ature of a star, computed from models of stellar atmo- 
spheres. For Cepheids, the values used are k = 0.064 
and a = 0.465 (|Abbottlil9"82l ). The luminosity, mass and 
radius are denoted L, M, and R respectively while a e 
is the electron scattering opacity, and Vth is the thermal 
velocity. The function T e = a e L / '(AircGM) is the ratio 
of the acceleration due to continuum driving and gravity, 
while the function £ is the sum of the acceleration due 
to pulsation and shocks, 
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The first term of Equation [2] represents the acceleration 
due to pulsation where AR is the amplitude of radius 
variation, to = 2n/P is the angular frequency of pulsation 
and t is the time. The second term of ( is the acceler- 
ation due to shocks, where P is the period of pulsation, 
dus/d(j) is the change of velocity with respect to phase 
of the shocks at the surface for a model of the atmo- 
sphere of S Cep ()Fokin et al.l Il996h . L$ : P$, and M$ are 
the luminosity, pulsation period and mass of the model 
of (5 Cep, and L and M is the luminosity and mass of 



Cepheid Mass Loss and Metallicity 



3 



the Cepheid. The function £ and shock acceleration are 
derived in Paper I. 

This approximate solution for the mass-loss rate of a 
pulsating star can be used to evaluate the dependence on 
metallicity. Consider first the dependence of mass loss 
on metallicity for winds driven by radiation alone. In 
this case the mass-loss rate is given by Equation Q] with 
£ = at all phases of pulsation. For Cepheids evolving 
on the second crossing of the instability strip the lumi- 
nosity is proportional to the mass , with the same depen- 
dence for different metallicities ([Castellani et al.l Il992t 
iGirardi et alJ 12000). The insta bility strip is som ewhat 
hotter for smaller metallicities (Bono et alJl2000f ). but 
there is significant overlap for the metallicities in ques- 
tion. Therefore the temperature range does not con- 
tribute significantly. If we consider two Cepheids with 
differing metallicities but same mass and luminosity then 
from Equation [1] with £ = the comparison of the mass- 
loss rates for two metallicities Z\ and Zi is 
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This implies that Cepheids in the LMC and SMC would 
have lower mass-loss rates then those in the Milky Way 
as the relative metallicity for the LMC is 0.4 and for the 
SMC is 0.2. Because a = 0.465, Mlmc/M mw = 0.139 
and Msmc/^mw = 0.031. Thus the radiative-driven 
mass-loss rates in the LMC and SMC are very small 
relative to the Milky Way. However the metallicity is 
important in the calculation of pulsation-driven mass- 
loss rates as well as radiative-driven mass-loss rates. If 
we again consider Equation [T] to describe the ratio of 
mass-loss rates for two different metallicities then 
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where the terms ashock and a p are the shock and the pul- 
sation acceleration terms of the function £ respectively, 
and g e g is the effective gravity, GM(1 — T e ). The ratio 
of the shock acceleration to the effective gravity is 
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and the ratio of the pulsation acceleration and the cffcc 
tive gravity is 
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The ratio of the analytic mass-loss rates can be used to 
test the behavior of a Magellanic Cloud Cepheid relative 
to a Galactic Cepheid, and to investigate if pulsation- 
driven mass loss is significant for lower metallicities. To 



do this, Equation [5] must be rewritten such that all vari- 
ables are in relative non-dimensional units. Consider two 
Cepheids, labeled 1 and 2, and express the ratio of the 
pulsation acceleration for the Cepheid number 2 be in 
terms of the effective gravity, called a p . Second, write 
the balance of force for Cepheid number 2 as the non- 
dimensional parameter F = 1 — a p — ashock,2/ffcff,2- Also 
it is assumed that F e ~ and can be ignored in the cal- 
culation because the electron number is small for stars 
with effective temperatures ranging from 4000 to 6000X. 
The small electron number causes an acceleration that is 
small relative to the gravity. The acceleration due to 
pulsation for Cepheid 1 relative to Cepheid 2 is 
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and the ratio of the shock acceleration and the effective 
gravity is 

-1.29 
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Substituting these results into Equation [5] and writing all 
quantities with a subscript of 1 in units of the quantities 
denoted with a subscript of 2, the relative mass-loss rate 
of Cepheid 1 is 
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Equation [TD] can be used to probe the mass-loss rates for 
low metallicity Cepheids as a function of relative period 
with the following free parameters: a p , F, Mi, Ri, and 
ARi. 

If the metallicity Z\ is less than unity then the relative 
mass-loss rates, M%, tend to be smaller based on the 
explicit dependence of metallicity. However, the ratio of 
mass-loss rates may be greater than or equal to unity 
if the term inside the curly brackets is significantly less 
than one. This term is significantly less than one only 
when 
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and furthermore this term must be smaller than the bal- 
ance of forces, F. There are many ways for this function 
to be small, one of which is the relative period. If Pi < 1 
and Ri and Mi are approximately unity, then this func- 
tion is less than one; however how much less than one 
Pi needs to be will depend on the values of a p and F. 
The value of a p , the ratio of the acceleration due to pul- 
sation to the effective gravity, is generally in the range 
to 0.1 for Galactic Cepheids. Therefore a p plays only a 
small role in affecting the mass-loss rate, implying that 
the shock acceleration term is the dominate term, which 
means the parameter F needs to be small to satisfy Equa- 
tion [TlJ A small value of F suggests that the balance of 
forces for the reference Galactic Cepheid is small and the 
mass-loss rate is large. 
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Fig. 1. — The mass-loss rate of a LMC Cepheid relative to a similar Cepheid with solar metallicity as a function of period and the effect 
of varying the other parameters. When the parameters are held constant, they are Ri = 1, Mi = 1, AR± = 1, F = 5%, and a p = 5%. 



The inequality given by Equation II II is satisfied when 
the relative radius is greater than unity or when the mass 
is less than unity, but this also requires that the value 
of F be small. Varying the amplitude of radius pulsa- 
tion affects the result as well; if the amplitude is greater 
than unity for lower metallicity then the pulsation terms 
increases linearly and thus decreases Equation [TTJ The 
result of varying the values of a p , F, R\, Mi, and Ai?i 
in Equation 1101 gives the mass-loss rate of Cepheid 1, 
Mi , shown in Figure Q] as a function of Pi for the rela- 
tive metallicity of the LMC to the Milky Way. From the 
plots, it is clear that the relative mass-loss rate can be- 
come very large even at lower metallicity. Applying the 
same analysis for the mass loss of SMC Cepheids relative 
to Galactic Cepheids the curves shows similar behavior. 
The effect of the metallicity in the analysis is to shift 
the curve up and down along the mass-loss rate axis. 
For a relative mass-loss rate of unity the relative period 
is approximately the same for both Cepheids with LMC 
metallicity and Cepheids with SMC metallicity since it 
is in the nonlinear regime given by the term in the curly 
brackets in Equation [T0l 

The mass-loss rates are larger at lower metallicity if 
the shocks are more efficient. These shocks are gener- 
ated in the hydrogen and helium partial ionization zones 
of Cepheids and at lower metallicity the partial ioniza- 
tion zones produce more energy when the layers become 
ionized. This is because of the larger fractional mass 
hydrogen and helium. Therefore the shocks have more 
momentum to input into the wind. The lower metallicity 
and possibly larger abundanc e of helium also affe ct the 
global properties of Cepheids (|Marconi et alJ l2005) , such 
as the pulsation period, and these metallicity-dependent 
global properties appear in the shock acceleration for- 
mula making it implicitly metallicity dependent. 



3. COMPARISON OF THE ANALYTIC PREDICTIONS 
WITH OBSERVATIONS 

It has been argued that pulsation-driven mass loss gen- 
erally decreases for lower metallicities, but there are com- 
binations of parameters that predict larger relative mass- 
loss rates. One might ask what the parameters describing 
Magellanic Cloud Cepheids would predict. 

The average relative period of the LMC and SMC 
Cepheids can be determined using the metallicity cor- 
rection. The metallicity correction is a constant that 
is added to the Period-Luminosity relation to account 
for the fact that Cepheids with different metallicities 
have different luminosities for the same pulsation pe- 
riod. The Hu bble Key Project on the Extragalactic 
Distance Scale (|Kennicutt et al.lll998f) used the Period- 
Luminosity relation to determine the distance to galax- 
ies in the local group, but it required a correction for 
the metallicity of approximately 8{m — M)/5[M/H] = 
—0.25 mag/dex. Thi s correction factor has been veri- 
fied b y other studies jT ammaiin et al. 2003; Saka i et all 
1 2004 IGroenewegen etaf] 12004: rRomaniello et al.l 120051 : 
iGroenewegen fc Salarisll2003l : lGieren et al.ll2005f ) but the 
exact value is uncertain. For instance, fecent the- 
oretical studies using non-linear hydrodynamic mod- 
els find the period-luminosity relation depends signif- 
icantly on the metallicity as well as on the helium 
abundance { Fiorcutin o et alJ 12002: Ma rconi et al.ll2005h . 
ISasselov et alJ (|1997t ) argued that the metallicity correc- 
tion, <5/i, could be written as a function of Z, 



Sfi = 0.44 log 



(—) 



(12) 



Combining this with the period-luminosity relation for 
the LMC, 



M v = -2.760(logP-l) -4.218, (13) 
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(|Freedman etalll200l . one finds the period ratio for 
the LMC and SMC relative to Milky Way Cepheids 
with the same luminosity to be Plmc /-Pmw ~ 1.16 and 
-Psmc/-Fmw ~ 1-30. 

The relative mass of Cepheids in the LMC and 
SMC can be estimated by comparing the mass dis- 
crepancy in the Mage llanic Clouds and the Milky Way 
(|Keller fc Woodl 120061) . The mass discrepancy is a func- 
tion of metallicity implying that for the same luminosity 
the mass of Cepheids is a function of metallicity. There- 
fore the mass of SMC Cepheids may be up to approxi- 
mately 3% smaller than LMC Cepheids, which in turn is 
about 3% smaller than Galactic Cepheids. 

The values of the relative radius and relative amplitude 
of radius variation are difficult to approximate. Radial 
velocity observations suggest that the radii of Cepheids 
in the Magellanic Clouds follo w a similar Period-Radiu s 
relation as Galactic Cepheids (|Storm et all |2004 [20051 ). 
Therefore varying these parameters can provide a test of 
mass loss with the period ratios quoted for the LMC and 
SMC. Figure [3 shows the relative mass-loss rate for LMC 
(Left Panel) and SMC (Right Panel) Cepheids as a func- 
tion of the mean radius for corresponding relative period. 
The analysis assumes the value of F is 0.01, ARi = 1.15, 
and a p — 0.05. Each curve represents the possible rela- 
tive mass range suggested by mass discrepancy studies, 
0.96 to 1 for the LMC and 0.92 to 1 for the SMC. The 
results are insensitive to small variations of F, AR%, and 
a p , implying that mass loss in the Clouds is significant if 
the radii of Cepheids are about 10% and 16% larger in 
the LMC and SMC respectively than Galactic Cepheids 
with the same luminosity. This may be tested by mea- 
suring t he angular diameters of Cepheids using interfer- 
ometry, iMourard et al.1 (|2008f l argued that it is possible 
to measure these angular diameters with the Very Large 
Telescope Interferometer using differential interferome- 
try techniques. 

The plausibility of this required larger radius may be 
checked by considering the observed Period-Radiu s (PR) 
relation ((Gieren et all 119991 : [G roencwegc rJ l200l . The 
PR relation for Galactic, LMC and SMC Cepheids is 

log R = 0.68 log P+ 1.146, (14) 

and for the period ratios of 1.16 and 1.3 the relative LMC 
and SMC Cepheid radius is 1.1 and 1.2 respectively. This 
implies that it is possible that LMC and SMC Cepheids 
have larger mass-loss rates than Galactic Cepheids when 
the mass loss of the Galactic Cepheid is significant. 

The metallicity correction that is determined using 
theoretical models suggests that metal-poor Cepheids 
are more luminous for the same period. In other words, 
the periods of m etal-poor Cepheids ar e shorter for the 
same luminosity (jFiorentino et al.ll2002h . In this case the 
value of Pi is less than unity. The relative mass-loss rate 
will be greater than unity for LMC and SMC Cepheids 
as long as the other parameters are i?i > 1, ARi > 1 
and Mi < 1. 

This analysis has shown that it is possible for metal- 
poor Cepheids to have significant mass loss, similar to 
and even greater than the mass-loss rates for Galactic 
Cepheids for both relative period regimes suggested by 
theoretical and observed values of the metallicity correc- 
tion of the Period-Luminosity relation. 



4. PREDICTING MASS-LOSS RATES FOR THEORETICAL 
MODEL CEPHEIDS 

The previous analysis is suggestive, but it can be 
carried further to determine, quantitatively, how much 
larger the mass-loss rates are for LMC and SMC 
Cepheids, at least in the regime where the periods of 
LMC and SMC Cepheids are smaller than those of Milky 
Way Cepheids for the same luminosity. To compute the 
mass-loss rate for a Cepheid, it is necessary to know 
the following parameters that describe the Cepheid: the 
mass, luminosity, radius, and pulsation period, as well 
as the amplitudes of the variation of the luminosity 
and radius. While these quantit ies have been deter- 
mine d for many Galactic Cepheids (Moskalik fc Goryriyal 
120051 ). this has not been done for LMC and SMC 
Cepheids. Therefo re theoretical models of Cepheids from 
iBono et all (|2000l ) are used to predict mass-loss rates. 
The method of solving the mass loss in the pulsating 
case is derived in Paper I. 

The models include all the necessary information to 
calculate pulsation-driven mass-loss rates using the en- 
hanced CAK method. There are two types of models 
describing the mass-luminosity relations: one is based 
on canonical stellar evolution models, and the second as- 
sumes convective core overshoot, which results in a larger 
luminosity for the same mass as the canonical models. 
This does not mean the models necessarily represent only 
convective core overshoot; for instance the steeper mass- 
luminosity relation could also represent mass loss during 
stages of evolution before the second crossing of the in- 
stability strip. The masses of the models are 5, 7, 9 and 
IIMq and span the temperature range of the instability 
strip. 

The mass-loss rates for model Cepheids in the Milky 
Way, LMC and SMC are shown in Figure [3] as a func- 
tion of pulsation period for the canonical models, along 
with the ratio of the pulsation mass-loss rates and the 
radiative-driven mass-loss rates. The ratio, the enhance- 
ment of mass loss, is a measure of the dependence of the 
wind on the pulsation plus shock terms in the mass-loss 
calculation. The points labeled as LMC or SMC Maxi- 
mum are based on model Cepheids where the pulsation 
period is so small that the sum of the acceleration due to 
pulsation and shocks is greater than the effective gravity 
of the Cepheid. In this limit the calculation of the mass- 
loss rate becomes imaginary. These points are calculated 
by using the same parameters from the model Cepheid 
but increasing the period by steps of 0.01 day until the 
calculation is stable. The physical justification for this 
is based on how the rate of change of period depends on 
mass loss. It was shown in Paper I that mass loss acts to 
increase the rate of period change, making it more posi- 
tive. For most Cepheids, the contribution to the period 
change due to mass loss is small compared to the contri- 
bution due to evolution, but for Cepheids near the blue 
edge of the instability strip the effect of evolution on pe- 
riod change decreases, making the contribution of mass 
loss more significant. The period is also smaller, further 
increasing the mass-loss rate and the effect of mass loss 
on the period change. 

The mass-loss rates for the SMC models shown in 
Figure [3] are larger than that for LMC and Milky Way 
at periods about 3.5 days, denoted by LMC and SMC 
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Fig. 2. — The relative mass-loss rates for LMC (Left Panel) and SMC (Right Panel) Cepheid as a function of relative radius with period 
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Fig. 3. — (Lef t Panel) The predicted mass-loss rates of the theoretical Cepheid models with a canonical mass-luminosity relation 
(Bono ct al. 2000). (Right Panel) The enhancement of the mass loss as given by the ratio of the pulsation-driven mass-loss rate to the 
radiative— driven mass— loss rate that ignores pulsation effects. The points given in the key as LMC/SMC Maximum are explained in the 
text. 

Maximum in the figures. This is also true at periods 
near 20 days. These periods correspond to masses of 5 
and 9Mq respectively. The increase of mass-loss rates 
with decreasing metallicity is clear in these two cases, 
but for the other models the effect is not as obvious. It 
is interesting that the long period Cepheids have lower 
mass-loss rates with decreasing metallicity. This is con- 
sistent with the observation that the mass discrepancy of 
Ceph eids decrease with i ncreasing period for the Milky 
Way (jCaputo et al.ll2005l ), and it suggests similar behav- 
ior for LMC and SMC Cepheids. The mass-loss enhance- 
ment, in the right panel of Figure [3j shows that the mass 
loss in the SMC Cepheids is more sensitive to the effects 
of pulsation for Cepheids in the LMC and Milky Way; for 
Cepheids with periods near 10 days the SMC mass loss is 
enhanced by a factor of 10 relative to the LMC. This ver- 
ifies the the argument that when mass loss is enhanced 
by pulsation for Galactic Cepheids, the mass loss is fur- 
ther amplified for similar Cepheids at the same position 
on the instability strip but with lower metallicities. 

The analysis is repeated with the steeper mass- 
luminosity relation associated with convective core over- 
shoot and is shown in Figure [U The dependence of 
pulsation-driven mass loss on metallicity is more clear 



in this case. The mass-loss rates for SMC Cepheids may 
be larger than those of the LMC and the Milky Way 
for periods up to about 20 days which is consistent with 
the result shown in Figure [3] For periods greater than 
20 days, the mass-loss rate decreases as a function of 
metallicity where pulsation and shocks are not as effec- 
tive. The mass-loss enhancement, shown in the right 
panel of Figure [4J implies pulsation is more efficient for 
driving mass loss as metallicity decreases. The mass-loss 
rates are also largely independent of the mass-luminosity 
relation for the theoretical Cepheid models. 

5. PREDICTIONS OF INFRARED EXCESS DUE TO MASS 

LOSS 

It has been hypothesized that mass loss is important 
in LMC and SMC Cepheids. In Paper I, the circum- 
stellar shells that have been observed as infrared excess 
in Galactic Cepheids were modeled as dust forming in 
the winds of Cepheids. This luminosity depends on the 
fraction of the total amount of gas that forms dust, the 
dust-to-gas ratio, which for the Milky Way is approxi- 
mately 1/100. This quantity is dependent on the compo- 
sition of the gas, meaning the dust-to-gas ratio is smaller 
for the lower metallicity of the LMC and SMC. For this 
work the dust-to-gas ratio is assumed to scale linearly 
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Fig. 4. — Same as Figure[3]but for Cepheid models assuming a convective overshoot mass-luminosity relation. 
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with the relative metallicity, meaning that the ratio for 
the LMC is 1/250 and SMC is 1/500. However these 
are m aximum values for t he du st-to-gas ratios; for in- 
stance IClavton k. Martini (|1985t ) found the LMC dust- 
to-gas ratio is about 1 /400 and iWeingartner k. Draind 
(2001) determined that the ratio is about 1/1000. These 
chosen dust-to-gas values are the maximum values but 
are uncertain by only about a factor of 5. 

If mass loss is important in LMC and SMC Cepheids 
then the winds should produce observable infrared ex- 
cess. A Color-Magnitude Diagram, My versus My — 
M 3 . 6nmi is shown in Figure [5] for the Cepheid models 
based on the canonical mass-luminosity relation. Along 
with that plot, the difference of the colors due to the star 
plus circumstellar shell relative to the star alone is also 
shown in the right panel. Because the luminosity of a 
circumstellar shell is due to dust, it will have minimal 
effect on the value My, but the dust contributes signifi- 
cantly to the 3.6 /xm luminosity. Therefore color excess 
increases if the mass-loss rate is large. 

The color-excess analysis is repeated for the convec- 
tive core overshoot models and is shown in Figure El In 
this case, the color excess is more pronounced though for 
only a few points. In both cases, the color excess due 
to circumstellar shells is significant. If mass loss is im- 
portant in LMC and SMC Cepheids then observations of 



a large population of Cepheids should be able test this 
model. 

Infrared excess may also affect the infrared Period- 
Luminosity relation. The Period-Luminosity rela- 
tion for the LMC wa s rec ently presented by b oth 
iNgeow fc Kanburi (|2008ft and IFreedman et al.1 (I2008D i n 
IRAC bands using SAGE data (jMeixner et alj 120061) . 
The PL rel ation was also deter mined in J, H and 
K-bands bv iPersson et al.1 (I2004D fo r a sample of 92 
Cephe ids. In lNgeow fc Kanb ur (2008) and lPersson et al.1 
(2004) blending is an important issue because of the 
lower resolution of t he infrared observations; while 
IFreedman et al.l (l2008f) avoids th e issue by using a sub- 
sample from lPersson et ail (|2004) . Blending acts to make 
the Cepheid appear more luminous in the infrared, just 
as mass loss would do, making the two difficult to dis- 
tinguish from photometry. In Figure [71 the luminosities 
of the theoretical Cepheids, including infrared excess due 
to mass loss are plotted as a function of period for the 
Galactic, LMC and SMC sam ple at 3.6 ixm. The 3.6 
fim PL relation for the LMC (jNgeow fc Kanburi [2001) 
is also plotted with dotted lines representing 3er devia- 
tion. For the model LMC Cepheids with the canonical 
mass-luminosity relation, there are two regimes where 
the luminosity is 3tr larger than the luminosity from the 
observed PL relation. There are two Cepheids at a pe- 
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riod of about 15 days and the two near four days having 
the maximum possible mass-loss rate with a luminosity 
four magnitudes larger than that given by the PL rela- 
tion. The models with convective overshooting, shown 
in the Right Panel of Figure \7\ show a similar behavior 
at 15 days, but at shorter periods the luminosity of the 
LMC Cepheids are consistent with the PL re lation. 

In the analysis of iNgeow fc Kanburl (|2008l ). a number 
of Cepheids are removed from the determination of the 
infrared PL relation because the IR brightness is too 
large. This is presumably due to blending, however in- 
frared excess due to mass loss has the same effect. The 
predicted luminosity of the model Cepheids is consistent 
with the range of luminosities of these outliers, providing 
possible observational agreement. 

The computation of mass-loss rates of the theoretical 
Cepheid models predicts that at some points on the sec- 
ond crossing of the instability strip the mass-loss rates 
are increased by up to three orders of magnitude. This 
significant increase agrees with the result in the section 
2. 

6. MASS LOSS AND THE MASS DISCREPANCY 

It has been shown that mass loss increases as metal- 
licity decreases for theoretical models of Cepheids. This 



result is consistent with the observatio ns of mass dis- 
crepancy as a fu nction of metallicity (|Keller fc Woodl 
120061 : [Keller 20Q8|). Furthermore the dependence of mass 
loss on metallicity is found to change for larger periods 
where the mass-loss enhancement becomes less signifi- 
cant. Along with the fact that more massive Cepheids 
spend less time on the instability strip than less massive 
Cepheids, this result is potentially consistent with the ob- 
serva tion that the mass discrepancy is constant (jKellerl 
l2008f) or decreases (|Caputo et al.ll2005h as a function of 
mass. The model of mass loss being driven by radial 
pulsation and shocks is consistent with the behavior of 
mass discrepancy between pulsation and evolution calcu- 
lations as a function of both mass and metallicity. But 
can mass loss predict the actual measured differences of 
mass? 

For mass loss to be the solution it must be able 
to account for a difference of approximately 1M@ for 
Cepheids with an evolutionary mass of 5Mq and - 
2Mq for Cepheids with an evolutionary mass 11M Q in 
the Milky Way. The lower limit represents the mass dis- 
cre pancy with th e depen dence on mass that is predicted 
bv ICaputo et all (|2005f ) and the upper l imit is due t o 
the constant mass discrepancy argued by Keller; (;2008;). 



Cepheid Mass Loss and Metallicity 



9 



For LMC Cepheids the mass difference is approximately 
1 - 1.5M© for 5M Cepheids and - 2.5M Q for 11M© 
Cepheids, while for SMC Cepheids the mass difference is 
1 - 2Mq and - 3M respectively. 

A timescale is required to determine the average mass- 
loss rates needed to account for mass discrepancy. This 
timescale is assumed to be the evolutionary lifetime on 
the second crossing: about 25 Myr for a 5M Galactic 
Cepheid and abou t 2 Myr for an 11M Galactic Cepheid 
(|Bono et al.ll2006f) . Furthermore, it is assumed that the 
evolutionary timescales are similar in the LMC and SMC. 
In that 5M Q Cepheid would need to lose mass at 

an average rate of 4 x 10~ 8 in the Galaxy, 6 x 10 -8 in the 
LMC and 8 x lO~ 8 M /yr in the SMC, to account for the 
maximum mass difference. For an 11 Mq Cepheid, the 
average mass-loss rate needs tobeO-lxlO -6 , 0-1. 2 x 
10~ 6 and - 1.5 x 10~ 6 Af Q /yr. According to the theory 
of radiative-driven mass loss, which dominates over the 
pulsation effects at long period, the mass-loss rate is of 
order 1O _8 M , far less than the maximum value of the 
average mass-loss rate that is required. The mass-loss 
rate for I Car has been estima ted to be < 2 x 10~ 8 M@/yr 
(jBohm-Vitense fc Lovelll994t ). and for RS Pup the rate 
is estimated to be < 3.5 x l(r 6 M Q /V (|Deasvlll988r ). 

The challenge for explaining the mass discrepancy with 
mass loss is thus the large-mass regime. However, this 
challenge disappears if the mass di screpancy is a decreas - 
ing function of mass, as found bv ICaputo etaH (|2005h . 
The mass discrepancy does not need to be zero at large 
mass, but it must be smaller than the 17% argued by 
iKelTin ([2008) to agree with predictions from pulsation- 
driven mass loss. If the mass discrepancy is smaller for 
large mass Cepheids, then an average mass-loss rate of 
the order lO _8 M /yr may resolve the difference. 

At smaller mass, M < 9M , the theory of pulsation- 
driven mass loss predicts mass-loss rates of order 
lO -8 M /?/r consistent with the average required mass- 
loss rate for the mass discrepancy. It is not obvious, 
however, that the average mass-loss rate for Cepheids 
over the first and second crossing is this large. The pre- 
dicte d mass-loss rates of the model Galactic Cepheids 
from iBono et al.l (|2000D is generally too small, but there 
are a number of examples of significant mass loss in short 
period Cepheids in Paper I. 

In the LMC and SMC the mass discrepancy at smaller 
mass is 17 - 25%, and the required average mass-loss rate 
to explain this is 6 - 8 x lCr 8 Af /?/r. At 10 - 20 day 
periods the mass-loss rates of the theoretical Cepheid 
models is greater or equal to this required average, as 
shown in Figure [3j but only for two of the points of on 
the evolutionary track of the LMC and SMC Cepheid. 
For the shorter period LMC/SMC Cepheids the mass 
loss was shown to be unstable for some of the models, 
but this was solved by using slightly larger values for the 
period of pulsation for a maximum value of the mass-loss 
rate. These maximum values are of the order 1O _7 M 
implying that it is possible to account for the required 
amount of mass loss. 

It is worth noting that pulsation-driven mass loss pre- 
dicts that the mass loss decreases with metallicity at long 
periods. If mass loss is the solution to the mass discrep- 
ancy then the mass discrepancy must decrease as a func- 
tion of mass at a larger rate for the LMC a nd SMC than 
for the Milky Way, against the assertion of lKellerl (|2008l ) 



that the mass discrepancy is constant with mass. The 
mass discrepancy may be constant for periods up to 20 
days where the mass loss is dominated by pulsation and 
shocks. At longer periods, the mass loss depends pri- 
marily on radiative driving, which decreases as (Z/Zq) 2 7 
thus implying a smaller mass discrepancy in the SMC 
and LMC compared to the Milky Way. 

The results of this work imply it is possible for mass 
loss to account for the mass discrepancy but the evi- 
dence is circumstantial. However, mass loss accounts 
for at least some of the mass discrepancy according to 
this theory. In that case, it could be argued the remain- 
ing mass discrepancy is due to convective core overshoot 
which in turn affects mass loss. The most significant ef- 
fect of convective core overshoot on the global parameters 
is the ratio of the luminosity and mass, where this ratio 
is larger for more overshoot. This larger ratio causes a 
larger contribution to the acceleration of the wind via the 
continuum opacity which in turn increases the mass-loss 
rate of all Cepheids. By comparing the mass-loss rates 
in Figures [3] and |U it is reasonable to argue that mass 
loss is important even if the helium core is larger due to 
convective core overshoot in Cepheid progenitors. The 
two possible explanations may be differentiated by deter- 
mining the mass discrepancy for long period (P > 30d), 
massive Cepheids; if the discrepancy is smaller at large 
masses then mass loss is a probable solution. 

7. CONCLUSIONS 

In this article, we have investigated the hypothesis that 
mass loss is important for the lower metallicity Cepheids 
in the LMC and SMC. This has been done by rewrit- 
ing the analytic method for computing pulsation-driven 
mass-loss rates as a ratio of mass-loss rates for two 
Cepheids that are alike in effective temperature and lu- 
minosity but with different metallicities, and allowing 
the other parameters to vary, such as the radius and the 
period of pulsation. This leads to a parameter space 
that can be tested to understand the potential behav- 
ior of Cepheids at lower metallicity. Over the major- 
ity of the parameter space the relative mass-loss of the 
two Cepheids is less than unity, reflecting the explicit 
dependence on metallicity. However, there exists two 
regimes where the relative mass-loss may be larger for 
lower metallicity. The first regime is where the relative 
mass is less than unity and the relative radius is larger 
than unity, and thus the ratio of the pulsation periods 
may be greater than unity. The other regime is where the 
pulsation period of the lower metallicity Cepheid is less 
than that of the Cepheid with solar metallicity and the 
ratios of the mass and radius are approximately unity. 
There is an additional parameter constraining this re- 
sult. The mass loss is only larger for the lower metallic- 
ity Cepheid if the mass loss is significant for the similar 
solar-metallicity Cepheid. In other words, if mass loss is 
significantly amplified by pulsation for the solar metallic- 
ity Cepheid then the mass loss may be further amplified 
by pulsation for the lower metallicity Cepheid. 

We tested if the parameter space that caused the pul- 
sation amplification of mass loss is consistent with with 
properties of LMC and SMC Cepheids. This was done 
using the metallicity correction of the Period-Luminosity 
relation to determine periods of LMC and SMC Cepheids 
relative to Galactic Cepheids with the same luminos- 



10 



Neilson & Lester 



ity. Also we used the Period-Radius relation to deter- 
mine how the radius changes for the lower metallicity 
Cepheids. The metallicity correction suggests that lower 
metallicity Cepheids have longer pulsations periods than 
Galactic Cepheids and the longer pulsation periods im- 
ply that the radii of lower metallicity Cepheids are larger 
than Galactic Cepheids for the same luminosity. These 
behaviors are consistent with mass-loss rates being fur- 
ther amplified implying mass loss is important for LMC 
and SMC Cepheids. 

We have calculated mass-loss rates for theoretical 
models of Large, and Small Magellanic Cloud, and 
Galactic Cepheids using the method presented in Pa- 
per I. It should be noted that the theoretical mod- 
els use stellar evolution calculations that do not in- 
clude mass loss in earlier stages of evolution. This af- 
fects the Cepheid mass-luminosity relation. It is not 
clear how the pulsation-driven mass-loss rates of mas- 
sive (M > WMq) Cepheids would cha nge if we us e 
the non-linear mass-luminosity relation of lKellerl (2008). 
The non-linear mass-luminosity relation suggests that at 
large mass the lumin osity is less fo r a giv en mass than 
what is predicted bv ICaputo et al.l ()2005l ). This would 
suggest that pulsation-driven mass-loss rate would be 
less, but it is unclear what the pulsation period would 
be in this case. The pulsation period of a Cepheid fol- 
lowing the non-linear mass-luminosity relation may be 
smaller than the period of a Cepheid following the linear 
mass-luminosity relation. If the period is smaller then 
the shocks that help drive mass-loss are more efficient 
and the mass-loss rate may increase. It would be very 
interesting to test this possibility. 

In the calculation of the mass-loss rates, it was found 
that the peak mass-loss rates in the second crossing of 
the instability strip in the SMC is almost a factor of 
10 larger than in the LMC, which in turn, is almost 10 
times larger than that in the Milky Way. Furthermore 
it was found that the period of about 10-20 days is the 
point where the dependence of mass loss on metallicity 
changes sign. This result may be interpreted as the mass 
loss being amplified significantly for the shorter period 
Cepheids and radiative driving being the dominant driv- 
ing mechanism in longer period Cepheids. The larger 
mass-loss rates for lower metallicity Cepheids is consis- 
tent with the conclusion that mass loss may be amplified 
more by pulsation at lower metallicity. 

It is believed that dust forms in the wind of a Cepheid, 
generating an infrared excess. From the predicted mass- 
loss rates of theoretical models of Cepheids, we compute 
Color Magnitude Diagrams as an observational predic- 
tion. The result, however, depends on both the dust-to- 
gas ratio and the mass-loss rate. The color excess com- 
puted here is thus the maximum value because the dust- 
to-gas ratio used in this work is the maximum value. 



The dust-to-gas ratio may be a few times to an order of 
magnitude less, implying the color excess may be about 
a factor of 2-3 less. The predicted 3.6 /im luminosity 
computed from the sum of the stellar and circumstel- 
lar shell luminosities was found to be consistent with 
the o bserved infrared Period-Luminosity relat i on at 3.6 
fim (|Ngeow fc Kanburl l2008t iFreedman et~aT1 120081) as 
well as possibly a g reeing with the set of Cepheids from 
iNgeow fc Kanburl (|2008f l that were considered outliers 
with much larger luminosities due to blending effects. 

The ability of mass loss to explain the mass discrep- 
ancy has also been explored. It is shown that pulsation- 
driven mass loss is a possible explanation, but there still 
is not enough information to say this with certainty. The 
most difficult aspect to explain is the 17% mass discrep- 
ancy at large mass where the predicted mass-loss rates 
are smaller than the necessary values of 10 -6 M@/yr. 
However, the 17% mass discrepancy is an upper limit 
and the mass discrepancy at large masses range from 
17% down to zero. At smaller masses the pulsation- 
driven mass loss is consistent with a mass discrepancy 
of 17%. For the LMC and SMC the mass-loss rates in- 
crease, implying mass loss can also explain the metallicity 
dependence of the mass discrepancy. 

The results of this work suggest that it is possible for 
mass-loss rates of lower metallicity Cepheids in the LMC 
and SMC to be greater than the mass-loss rates of sim- 
ilar Galactic Cepheids. This result is counterintuitive, 
one might expect mass loss to scale as the metallicity. 
However, the pulsation period, radius and other param- 
eters that describe a Cepheid and its mass loss also have 
dependencies on the metallicity, and the combination of 
these parameters and the dependencies cause the mass- 
loss rates to be potentially more significant. Further- 
more, the amount of mass loss predicted for Magellanic 
Cloud Cepheids is consistent with the observed behav- 
ior of the mass discrepancy of Cepheid s as a function of 
metallicity from lKeller fc Wood! (120061) as well as the de- 
pendence on mass lCaputo et all (|2005t ) . Therefore mass 
loss cannot be discounted as a possible contribution to 
the mass discrepancy. The amount of mass loss calcu- 
lated using the pulsation mass-loss model is significant 
but there is not enough information to argue mass loss 
alone resolves the mass discrepancy. This needs to be 
tested by with a large number of nonlinear pulsation 
models of Cepheids on both the first and second crossing 
to account for the effect of mass loss of the evolution of 
the period. 



The authors appreciate the helpful comments from the 
referee. HRN is grateful for funding from the Walter 
John Helm OGSST and the Walter C. Sumner Memorial 
Fellowship. 



REFERENCES 



Abbott, D. C. 1982, ApJ, 259, 282 

Beaulieu, J. P., Buchler, J. R., & Kollath, Z. 2001, A&A, 373, 164 
Bohm-Vitense, E., Evans, N. R., Carpenter, K., Albrow, M. D., 

Cottrell, P. L., Robinson, R., & Beck-Winchatz, B. 1998, ApJ, 

505, 903 

Bohm-Vitense, E. & Love, S. G. 1994, ApJ, 420, 401 

Bono, C, Caputo, F., &; Castellani, V. 2006, Memorie della Societa 

Astronomica Italiana, 77, 207 
Bono, C, Castellani, V., & Marconi, M. 2000, ApJ, 529, 293 



Brocato, E., Caputo, F., Castellani, V., Marconi, M., Sz Musella, 

I. 2004, AJ, 128, 1597 
Brunish, W. M. & Willson, L. A. 1989, in IAU Colloq. Ill: The 

Use of pulsating stars in fundamental problems of astronomy, ed. 

E. G. Schmidt, 252-+ 
Caputo, F., Bono, G., Fiorentino, G., Marconi, M., & Musella, I. 

2005, ApJ, 629, 1021 
Castellani, V., Chieffi, A., & Straniero, O. 1992, ApJS, 78, 517 
Castor, J. I., Abbott, D. C, & Klein, R. I. 1975, ApJ, 195, 157 



Ccpheid Mass Loss and Metallicity 



11 



Chiosi, E. & Vallenari, A. 2007, A&A, 466, 165 

Clayton, G. C. & Martin, P. G. 1985, ApJ, 288, 558 

Cordicr, D., Lebreton, Y., Goupil, M.-J., Lejeune, T., Beaulieu, 

J.-P., & Arenou, F. 2002, A&A, 392, 169 
Cox, A. N. 1980, ARA&A, 18, 15 

de Jager, C, Nieuwenhuijzen, H., & van der Hucht, K. A. 1988, 

A&AS, 72, 259 
Deasy, H. P. 1988, MNRAS, 231, 673 

Evans, N. R., Massa, D., Fullcrton, A., Sonneborn, G., & Iping, R. 

2006, ApJ, 647, 1387 
Fiorentino, G., Caputo, F., Marconi, M., & Musella, I. 2002, ApJ, 

576, 402 

Fokin, A. B., Gillet, D., & Breitfellner, M. G. 1996, A&A, 307, 503 
Frcedman, W. L., Madore, B. F., Gibson, B. K., Ferrarese, L., 
Kelson, D. D., Sakai, S., Mould, J. R., Kennicutt, Jr., R. C, 
Ford, H. C, Graham, J. A., Huchra, J. P., Hughes, S. M. G., 
Illingworth, G. D., Macri, L. M., & Stetson, P. B. 2001, ApJ, 
553, 47 

Freedman, W. L., Madore, B. F., Rigby, J., Persson, S. E., & 

Sturch, L. 2008, ApJ, 679, 71 
Gieren, W., Storm, J., Barnes, III, T. G., Fouque, P., Pietrzyhski, 

G., & Kienzle, F. 2005, ApJ, 627, 224 
Gieren, W. P., MofTett, T. J., & Barnes, III, T. G. 1999, ApJ, 512, 

553 

Girardi, L., Bressan, A., Bcrtelli, G., & Chiosi, C. 2000, A&AS, 
141, 371 

Groenewegen, M. A. T. 2007, A&A, 474, 975 

Groenewegen, M. A. T., Romanicllo, M., Primas, F., & Mottini, 

M. 2004, A&A, 420, 655 
Groenewegen, M. A. T. & Salaris, M. 2003, A&A, 410, 887 
Havlen, R. J. 1972, A&A, 16, 252 
Huang, R. Q. & Weigert, A. 1983, A&A, 127, 309 
Keller, S. C. 2008, ApJ, 677, 483 
Keller, S. C. & Wood, P. R. 2002, ApJ, 578, 144 
— . 2006, ApJ, 642, 834 

Kennicutt, Jr., R. C, Stetson, P. B., Saha, A., Kelson, D., Rawson, 
D. M., Sakai, S., Madore, B. F., Mould, J. R., Freedman, W. L., 
Bresolin, F., Ferrarese, L., Ford, H., Gibson, B. K., Graham, 
J. A., Han, M., Harding, P., Hoessel, J. G., Huchra, J. P., 
Hughes, S. M. G., Illingworth, G. D., Macri, L. M., Phelps, R. L., 
Silbcrmann, N. A., Turner, A. M., & Wood, P. R. 1998, ApJ, 498, 
181 

Kervella, P., Merand, A., Pcrrin, G., & Coude Du Foresto, V. 2006, 
A&A, 448, 623 

Leavitt, H. S. & Pickering, E. C. 1912, Harvard College 

Observatory Circular, 173, 1 
Marconi, M., Musella, I., & Fiorentino, G. 2005, ApJ, 632, 590 
McAlary, C. W. & Welch, D. L. 1986, AJ, 91, 1209 
Meixner, M., Gordon, K. D., Indebetouw, R., Hora, J. L., Whitney, 

B. , Blum, R., Reach, W., Bernard, J.-P., Meade, M., Babler, B., 
Engelbracht, C. W., For, B.-Q., Misselt, K., Vijh, U., Leitherer, 

C, Cohen, M., Churchwell, E. B., Boulanger, F., Frogel, J. A., 
Fukui, Y., Gallagher, J., Gorjian, V., Harris, J., Kelly, D., 
Kawamura, A., Kim, S., Latter, W. B., Madden, S., Markwick- 
Kcmper, C, Mizuno, A., Mizuno, N., Mould, J., Nota, A., Oey, 
M. S., Olsen, K., Onishi, T., Paladini, R., Panagia, N., Perez- 
Gonzalez, P., Shibai, H., Sato, S., Smith, L., Staveley-Smith, L., 
Tielens, A. G. G. M., Ueta, T., Dyk, S. V., Volk, K., Werner, 
M., & Zaritsky, D. 2006, AJ, 132, 2268 



Merand, A., Aufdenberg, J. P., Kervella, P., Foresto, V. C. d., ten 
Brummclaar, T. A., McAlistcr, H. A., Sturmann, L., Sturmann, 
J., & Turner, N. H. 2007, ApJ, 664, 1093 
Merand, A., Kervella, P., Coude Du Foresto, V., Pcrrin, G., 
Ridgway, S. T., Aufdenberg, J. P., Ten Brummclaar, T. A., 
McAlister, H. A., Sturmann, L., Sturmann, J., Turner, N. H., 
& Berger, D. H. 2006, A&A, 453, 155 
Moskalik, P., Buchlcr, J. R., & Marom, A. 1992, ApJ, 385, 685 
Moskalik, P. & Gorynya, N. A. 2005, Acta Astronomica, 55, 247 
Mourard, D., Nardctto, N., Lagarde, S., Petrov, R., Bonneau, D., 
& Millour, F. 2008, in The Power of Optical/IR Interferometry: 
Recent Scientific Results and 2nd Generation, ed. A. Richichi, 

F. Dclplanckc, F. Parcsce, & A. Chelli, 477-+ 

Mucciarelli, A., Ferraro, F. R., Origlia, L., & Fusi Pccci, F. 2007, 

AJ, 133, 2053 
Neilson, H. R. & Lester, J. B. 2008, ApJ, 684, 569 
Ngeow, C. & Kanbur, S. M. 2008, ApJ, 679, 76 
Persson, S. E., Madore, B. F., Krzemihski, W., Freedman, W. L., 

Roth, M., & Murphy, D. C. 2004, AJ, 128, 2239 
Ribas, I., Jordi, C, & Gimenez, A. 2000, MNRAS, 318, L55 
Rodrigues, L. L. & Bohm-Vitense, E. 1992, ApJ, 401, 695 
Rogers, F. J. & Iglesias, C. A. 1992, ApJS, 79, 507 
Romanicllo, M., Primas, F., Mottini, M., Groenewegen, M., Bono, 

G. , & Frangois, P. 2005, A&A, 429, L37 

Rosvick, J. M. & Vandcnberg, D. A. 1998, AJ, 115, 1516 
Sakai, S., Ferrarese, L., Kennicutt, Jr., R. C, & Saha, A. 2004, 
ApJ, 608, 42 

Sasselov, D. D., Beaulieu, J. P., Renault, C, Grison, P., Ferlet, R., 
Vidal-Madjar, A., Maurice, E., Prevot, L., Aubourg, E., Bareyre, 
P., Brehin, S., Coutures, C, Dclabrouillc, N., de Kat, J., Gros, 
M., Laurent, B., Lachieze-Rey, M., Lesquoy, E., Magneville, C, 
Milsztajn, A., Moscoso, L., Qucinncc, F., Rich, J., Spiro, M., 
Vigroux, L., Zylberajch, S., Ansari, R., Cavalier, F., Moniez, M., 
Gry, C, Guibert, J., Moreau, O., & Tajhmady, F. 1997, A&A, 
324, 471 

Schroder, K.-P., Pols, O. R., & Eggleton, P. P. 1997, MNRAS, 285, 
696 

Sebo, K. M. & Wood, P. R. 1995, ApJ, 449, 164 

Storm, J., Carney, B. W., Gieren, W. P., Fouque, P., Freedman, 

W. L., Madore, B. F., & Habgood, M. J. 2004, A&A, 415, 521 
Storm, J., Gieren, W. P., Fouque, P., Barnes, III, T. G., & Gomez, 

M. 2005, A&A, 440, 487 
Tammann, G. A., Sandage, A., & Reindl, B. 2003, A&A, 404, 423 
Testa, V., Marconi, M., Musella, I., Ripepi, V., Dall'Ora, M., 

Ferraro, F. R., Mucciarelli, A., Mateo, M., & Cote, P. 2007, A&A, 

462, 599 

Turner, D. G. & Evans, N. R. 1984, ApJ, 283, 254 
Weingartner, J. C. & Draine, B. T. 2001, ApJ, 548, 296 
Welch, D. L. & Duric, N. 1988, AJ, 95, 1794 

Willson, L. A. & Bowcn, G. H. 1986, in Lecture Notes in Physics, 
Berlin Springer Verlag, Vol. 254, Cool Stars, Stellar Systems and 
the Sun, ed. M. Zeilik & D. M. Gibson, 385-+ 



